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Abstract 

(N ; 

\ We suggest the possibility that the anomalies observed in the LSND experiment 

and the Gallium radioactive source experiments may be due to neutrino oscillations 
^ . generated by a large squared-mass difference of about 20 — 30 eV 2 . We consider 

<^> \ the simplest 3+1 four-neutrino scheme that can accommodate also the observed 

solar and atmospheric neutrino oscillations. We show that, in this framework, the 

disappearance of v e and in short-baseline neutrino oscillation experiments is 
■ mainly due to active-sterile transitions. The implications of the first MiniBooNE 

Qh! results, appeared after the completion of this paper, are discussed in an addendum. 

Neutrino oscillation experiments have shown that neutrinos are massive and mixed 
particles (see the reviews in Refs. [H El EJ, SI O E], [TJ, [8] ) . The observation of u e — > u^ T 
.£h ! oscillations with a squared-mass difference 

S; Ami OL ^8xlO- 5 eV 2 (1) 

in solar and reactor neutrino experiments and the observation of — > v T oscillations 
with a squared-mass difference 

Am 2 ATM ~ 2.5 x 10" 3 eV 2 (2) 

in atmospheric and accelerator neutrino experiments can be accommodated in the mini- 
mal framework of three- neutrino mixing, in which the three active flavor neutrinos u e , u^, 
and v T are superpositions of three massive neutrinos v\, u 2 , and z/ 3 . This three-neutrino 
mixing framework cannot explain through neutrino oscillations the LSND — > v e signal 
[HI EDI EED, H2] , which requires a squared-mass difference 

Am 2 SND > KT 1 eV 2 . (3) 
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Another anomaly observed in neutrino experiments is the disappearance of u e : s in the 
Gallium radioactive source experiments GALLEX [T3l [14"] and SAGE [T51 fT6l [T7]. These 
experiments are tests of solar neutrino detectors in which intense artificial 51 Cr and 37 Ar 
neutrino sources were placed near or inside the detectors. Both 51 Cr and 37 Ar decay 
through electron capture (e~ + 51 Cr — > 51 V + v e and e~ + 37 Ar — ■> 37 CI + v e ). The energies 
of the emitted neutrinos are, respectively E = 752.73 ± 0.24 keV and E = 813.5 ± 0.3 keV 
[T8] . The neutrinos emitted by the artificial sources were detected through the same 
reaction used for the detection of solar neutrinos [T9j : 

z/ e + 71 Ga^ 71 Ge + e-, (4) 

which has the low neutrino energy threshold E t h = 0.233 MeV. The weighted average 
value of the ratio R of measured and predicted 71 Ge production rates is [17] 

R = 0.88 ± 0.05 . (5) 

In Ref. [TTJ it has been suggested that this anomaly may be due to an overestimate of the 
theoretical cross section of the Gallium detection process in Eq. (T4]). However, a Gallium 
cross section rescaled by the factor in Eq. (0) leads to a significant deterioration of the 
fit of solar neutrino data [2D] . 

In this paper we consider the possibility that the anomaly observed in Gallium ra- 
dioactive source experiments is due to neutrino oscillations^. 

Since the neutrino path in the Gallium radioactive source experiments was of the 
order of 10 cm, an explanation of the observed disappearance of f e 's through neutrino os- 
cillations requires a large squared- mass difference Amg a . In fact, requiring an oscillation 
length Lq£. = 47rE/\AmQ a \ smaller than about 10 cm, we obtain 

Am^ a >20eV 2 . (6) 

Assuming CPT invariance, the survival probability of neutrinos and antineutrinos are 
equal. It follows that the disappearance of electron neutrinos at the level indicated by 
Gallium radioactive source experiments appears to be in contradiction with the results 
of reactor neutrino oscillation experiments (see the review in Ref. [22J), which did not 
observe any disappearance of electron antineutrinos with an average energy of about 4 
MeV at distances between about 10 and 100 m from the reactor source. Let us notice, 
however, that the oscillation length of reactor neutrinos implied by Eq. ([6]) is much shorter 
than 10 m: 

^actors < 4Qcm _ (?) 

Hence, in reactor neutrino experiments the oscillations due to Am^ are seen as an energy- 
independent suppression of the electron antineutrino flux by the factor in Eq. ©. A 
measurement of such a suppression requires a precise calculation of the absolute electron 
antineutrino flux produced in a reactoJE Since this calculation is rather difficult, it 

1 The results of the hrst GALLEX artificial 51 Cr source experiment [13] has been used in Ref. [21] in 
order to constrain the neutrino mixing parameters. 

2 Information on u e disappearance which is independent of the absolute flux calculation can be ob- 
tained through the measurement of the energy spectrum (assuming it to be known with small uncertain- 
ties) or the comparison between rates measured with different source-detector distances. In these cases, 
reactor neutrino experiments are not sensitive to oscillations generated by a squarcd-mass difference 
Am 2 > 2eV 2 , as one can see, for example, from Fig. 13a of Ref. [23"] . 
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is possible that its systematic uncertainties have been underestimated. Therefore, a v e 
disappearance in reactor neutrino oscillation experiments at the level indicated by Eq. (jHJ) 
with the oscillation length in Eq. (JTj) is not excluded with absolute certainty. 

In this paper, we consider the possibility that both the LSND and Gallium anomalies 
are due to neutrino oscillations, through the same large squared-mass difference 

A^LSND + Ga>20eV 2 . (8) 

We consider, for simplicity, a four-neutrino mixing scheme, in which the three active 
flavor neutrinos u e , u T , and one sterile neutrino v s are superpositions of four mas- 
sive neutrinos v\, z/ 2 , z/ 3 , and z/ 4 . This is the simplest scheme in which there are three 
independent squared-mass differences which can accommodate the hierarchy 

Am LSND+Ga • (9) 

Four-neutrino mixing have already been considered in many papers as the explanation 
of the LSND anomaly (see the reviews in Refs. [3], SI El E] ) - Here, we further constrain 
the allowed values of the large squared-mass difference and the mixing of the electron 
neutrino by requiring that AmL SND+Ga is responsible of both the LSND and Gallium 
anomalies. 

Since the so-called 2+2 schemes are disfavored by the data PlESlESlEZlEElESlEnil^, 
we consider a 3+1 scheme, in which there is a group of three neutrino masses which is 
separated from an isolated neutrino mass by the LSND + Ga mass splitting. In this case, 
we have 

Am SOL = Am 21> Am ATM = l Am 3ll - l Am 32l ' ( 10 ) 

Am LSND+Ga = l Am 4ll - l A "4l - l Am 4 3 l , i 11 ) 

where Am^ = m 2 , — m 2 . Furthermore, we take into account the upper limit 

m /3 <2.3eV (95% CL) , (12) 

obtained in the Mainz [31] and Troitzk [32] tritium /3-decay experiments on the effective 
electron neutrino mass [331 EH ES] 

4 

ml = Y,\Uek\ 2 m 2 k . (13) 
fc=i 

Since the three active flavor neutrinos must have large mixings with V\, i/ 2 , and u 3 in 
order to accommodate the observed oscillations due to Amg OL and Am\ TM , the only 
scheme allowed is the one in which U\, z/ 2 , and z/ 3 are light, with masses 

mi ,m 2 ,m 3 < 2.3 eV, (14) 

and z/4 is heavy, with mass 



m 4 ~ ^ Am2 gND+Ga > 4.5 eV . (15) 

In four-neutrino schemes, the average v e survival probability in the Gallium experi- 
ments is given by 

(P^ e > = l-isin 2 2tf Ga , (16) 
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where $Ga is an effective mixing angle. Interpreting R in Eq. (jSJ) as (P Ue -+ Ue ), we obtain 

sin 2 2$ Ga = 0.24 ±0.10. (17) 

In the 3+1 mixing schemes (see the review in Ref. [3]), the survival and transition prob- 
abilities in short-baseline neutrino oscillation experiments have the two-neutrino mixing 
forms (for a, (3 — e, /i, r, s) 

P Va ->u a = P, a ^ a = 1 - sin 2 2§ aa sin 2 (^§^) , (18) 



= Pu p ^u a = P^, = iV* a = sin 2 2d aP sin 2 ( ) (a ? 0) , (19) 



Am\ x L 

1e 

where L is the source-detector distance and the effective mixing angles are given by 

sin 2 2$ aa = 4 \U a4 \ 2 (1 - |f/ a4 | 2 ) , (20) 

sin 2 2<d a/3 = sin 2 2^ Q = 4 \U a4 \ 2 \U^\ 2 {a ^ 0) . (21) 

Therefore, we have 

sin 2 2tf Ga = sin 2 2$ ee = 4 \U eA \ 2 (l - \U ei \ 2 ) . (22) 

Taking into account that |[/e 4 | 2 is small, in order to accommodate the observed oscillations 
due to Amg OL and Am\ TM , we obtain, from Eqs. ffTTl) and f[2"2"j) . 

|£/ e4 | 2 ~ 1 sin 2 2tf Ga ~ 0.06 ±0.03. (23) 



In spite of the relatively heavy mass of z/j in Eq. (JToj) . the mixing of z/ e with 1/4 is not 
a problem for the bound in Eq. (Tl2l) on the effective electron neutrino mass in /5-decay 
experiments. In fact, the contribution of z/ 4 to is 

m /3 (z/ 4 ) = |^4|m 4 ~l.l±0.3eV( I ^) . (24) 

Therefore, the bound in Eq. ffT2l implies 

m 4 <10eV. (25) 

Taking into account also Eq. ([H]), we obtain the allowed range 

20 eV 2 < Am 2 SND+Ga < 100 eV 2 . (26) 

Let us now consider the LSND — > u e signal, which has been observed with the 
probability [12J 

PLSND = ( 2M ± Q 67 ± QA5 ^ x 1Q -3 _ 

Since we are considering large values of Am 2 jSND+Ga in the interval in Eq. f[2"B"]) . the 
transition probability measured in the LSND experiment is the averaged probability 

(iW e > = ^sin 2 2tf LSN D, (28) 



with the effective mixing angle given by (see Eq. (|2Tjl ) 



sin 2 2tf lsnd = sin 2 2^ = 4 | U e4 | 2 | | 2 . (29) 
Thus, from Eq. ( 1271) . we obtain 

sin 2 2$ LSND ~ (5.3 ± 1.6) x 1(T 3 . (30) 

Short-baseline ( v\ —*"v e oscillations generated by a large squared-mass difference have 
been recently searched, with negative results, in the CCFR [3B], KARMEN [57) . NuTeV 
[3"8"] , and NOMAD [3U] experiments. From Fig. 8 of Ref . [33] , one can see that the range 
of sin 2 2$lsnd in Eq. (ED]) is compatible with the results of the CCFR, KARMEN, and 
NuTeV experiments if the allowed interval of Am 2 SND+Ga in Eq. (f26l) is restricted to 

20eV 2 < Am 2 SND+Ga <30eV 2 . (31) 

In fact, although a combined analysis of all the relevant neutrino oscillations data yields a 
poor goodness of frtfl [HI S21 SSI SS SSI HE EH E] , if the fit is accepted, there is an allowed 
region in the sin 2 24 SND -Am^ SND+Ga at sin 2 2$ LSN d ~ (2 - 5) x 1(T 3 and Am^ SND+Ga 
in the interval in Eq. ( |3~TT) (see Fig. 3 of Ref.[46j and Fig. 4 of the first arXiv version 
of Ref. [40] ) . This allowed region appears to be in contradiction only with the exclusion 
curve obtained in the NOMAD experiment (see Fig. 8 of Ref. [39]). 

From now on, we consider Am 2 jSND+Ga in the interval in Eq. f[3"Tj) . 

The range of sin 2 2'(?LSND m Eq. fl3"Ul) and the determination of |?7 e4 | 2 in Eq. (|2"3"1) from 
the Gallium anomaly allow us to determine the allowed range of |£7^ 4 | 2 : from Eq. (12*9*1) . 
we obtain 

l9 sin 2 2?9TSND sin 2 snd , , 

U ^ = A\TT \2 ~ ■ 2o a ~ °- 02 ± °' 01 • 32 

4|(7 e4 | 2 sm z 2V Ga 

This small value of |f/ M 4| 2 implies that the effective mixing angle in short-baseline 
disappearance experiments is given by 

sin 2 2^ ~4|[/ m4 | 2 ~ 0.08 ±0.04. (33) 

This value of sin 2 2$ wt is compatible with the exclusion curves of the CDHSW [47] and 
CCFR [IB] — > oscillation experiments for Am 2 SND+Ga in the interval in Eq. (|3"TT) . 
It is interesting to notice that the results of the CDHSW disappearance experiment 
favor a Am 2 jSND+Ga in the range in Eq. (13"T1) . as remarked at the end of the appendix of 
Ref. EO]. 

Let us now consider the experimental bounds on z/ M — > v T and v e — > v T transitions 
obtained in short-baseline experiments (CHORUS [49], NOMAD [39] and CCFR [50] EI]). 
From Fig. 14 of Ref. [22], one can see that, for Am 2 SND+Ga in the range in Eq. (I3T1) . the 
effective mixing angles 

sin 2 2$ er = 4|f7 e4 | 2 |f/ T4 | 2 , sin 2 2^ T = 4 \U^\ 2 \U r4 \ 2 (34) 



3 The fit can be improved by introducing a second sterile neutrino [40] . in a so-called 3+2 mixing 
scheme. However, it seems to us that it is highly unlikely that the two large squared-mass differences 
happen to have just the right values in the small regions which are not excluded by the neutrino oscillation 
data. 
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are bounded by 

sin 2 2tf er < 1 x 1(T 1 , sin 2 2^ MT < 2 x 10~ 3 . (35) 

Taking into account the allowed ranges of \U ei \ 2 and \U^\ 2 in Eqs. ( 1231) and (1321) . the 
limit on sin 2 2$ eT does not give a significant bound, whereas the limit on sin 2 2$ Mr yields 



sin 2 2?V ^ sin 2 2?V sin 2 2$ Ga 
4|^ 4 | 2 ~ 4 sin 2 2tf LSND 



Ur±\ 2 = -^rPif * 7o Ga < 0.05 . (36) 



Therefore, also |t/ r 4| 2 is constrained to be small. It follows that 

\U s4 \ 2 = 1 - (\U e4 \ 2 + |[/ M | 2 + \U T4 \ 2 ) > 0.8 , (37) 

and the u e disappearance indicated by Gallium radioactive source experiments is mainly 
due to v e — ► z/ s transitions with an effective mixing angle given by 

sin 2 1d es = 4 \U e4 \ 2 \U S ±\ 2 ~ 0.2 ± 0.1 . (38) 

These transitions are compatible with the CCFR bound on v e — > z/ s transitions (Fig. 4 of 
Ref. [51]) for the effective squared mass difference AmL SND+Ga confined in the range in 
Eq. j3!}. 

The z/ e — > v s transitions due to Am 2 SND+Ga affect also solar neutrino experiments. 
Since the mixing of v s with u±, z/ 2 , and z/ 3 is small, in practice solar neutrino experiments 
should observe an average probability of disappearance of electron neutrinos into sterile 
neutrinos of the same value as the ratio R in Eq. (jSJ) measured the Gallium radioactive 
source experiments: 

(P Ue ^ s ) ~ X - sin 2 2$ es ~ 0.10 ± 0.05 . (39) 

It is interesting to notice that a comparison of the SNO Neutral-Current (NC) data with 
the Standard Solar Model (SSM) prediction is compatible with v e — > v s transitions at 
the level indicated in Eq. ( 1391) . although no evidence can be claimed, because of the large 
theoretical uncertainty of the SSM prediction. In fact, the equivalent flux of 8 B electron 
neutrinos measured in SNC0 through the NC reaction u + d — > p + n + v, which is equally 
sensitive to u e , v^, and u T , is 



$sno = j 4 94 ± o.21+°; 38 ) x 10 6 cm" 2 s" 1 . (40) 
This value can be compared with the BS05(GS98) [51] and TC04 [3S] SSM values 

$?| 05 = (5.69ig;||) x 10 6 cm" 2 s" 1 , (41) 

$JC04 = ( 5 31 ± q ^ x 1Q 6 cm -2 s -l ^ ^ 



leading to 



$SNO ( 

(^ e -«/ a )sNO+BS05 = 1 - .BS05 = °- 13 -(U7 > ( 43 ) 



^ 8 B 



4 The flux in Eq. (TH)]) has been measured in the phase II of the SNO experiment (also called "salt 
phase" ) , in which about 2 tons of NaCl have been added to the heavy water in order to improve the 
efficiency and precision of the NC measurement [53] . 
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$SNO 

(^.)SNO+TC04 = 1 " = °- 07 ± °- 13 • ( 44 ) 

One can see that, although the uncertainties are large, the tendency of the ratios in 
Eqs. (143 p and (144p is towards an agreement with the average probability of v e — > v s 
transitions in Eq. ( 1311 . 

The disappearance of v e due to v e — > v s transitions could affect the search for — > u e 
transitions in the MiniBooNE^] experiment [Hoi [oTf EB] , which has been designed to check 
the LSND anomaly. This is due to the fact that the MiniBooNE beam has a natural v e 
contamination of about 5 x 10~ 3 . Since the MiniBooNE detector is located at a distance 
of 541 m from the target and the energy spectrum of the beam ranges from about 
0.2 GeV to about 3 GeV, with a peak at about 0.6 GeV, it is convenient to write the 
oscillation length due to AmL SND+Ga as 

L ~ 120 m ( — \ f Am LSND+Ga \ ^ 

osc ~ U \GeVj{ 20 eV 2 J ' 1 j 

Hence, a AmL SND+Ga in the range in Eq. (13~T|) implies that the oscillation length is much 
shorter than the MiniBooNE source-detector distance and the flavor transitions are prac- 
tically constant over the energy spectrum. The effect on the u e spectrum at the detector 
is the superposition of two opposite and competitive contributions: a v e disappearance 
due to v e — ► v a oscillations with a relatively large mixing (see Eq. flHHj) ) and a v e ap- 
pearance due to — > v e with a relative small mixing (see Eq. f[3"01 ). Since the natural 
contamination of v e in the beam is at the percent level, the two opposite effects on 
the v e spectrum are competitive. 

The hypothesis of — > u e transitions driven by AmLSND+Ga ma y soon be tested 
at the T2K beam line (starting from 2009) with the near off-axis detector located at a 
distance of 280 m from the target. The neutrino energy in T2K is about the same as in 
MiniBooNE. With a systematic error on the electron neutrino flux cr(syst) ~ 5% the 90% 
CL. sensitivity to sin 2 2^ e is about 3 x 10~ 3 . 

The scenario under consideration implies also short-baseline — > v s oscillations gen- 
erated by Am 2 SND+Ga with the effective mixing angle 

sin 2 2$ MS = 4 \\J^f |£/" s4 | 2 ~ 0.08 ± 0.04 . (46) 

Since this value of sin 2 2$^ practically coincides with the value of sin 2 2$^ in Eq. fl3"3"j) 
and is much larger than the values of sin 2 2$^ e = sin 2 2 i # LSND in Eq. flHOl and sin 2 2$^ in 

Eq. (|33|) . the — > v s channel is the dominant cause of short-baseline ( z/^ disappearance. 

Optimal future experiments which could observe the large disappearance of ( z/g and 

due to active-sterile transitions and the ( z/^ { u e transitions due to AmL SND+Ga are: 
Beta-Beam experiments jSH] with a pure v e beam from nuclear decay (see the reviews in 
Refs. [HOI EI]); Neutrino Factory experiments with a beam composed of v e and z/ M , from 
/i + decay, or v e and from \i~ decay (see the review in Ref. [62J [60]); experiments with a 
v e beam produced in recoiless nuclear decay and detected in recoiless nuclear antineutrino 
capture [63J. 

5 The implications of the first MiniBooNE results, appeared after the completion of this paper, are 
discussed in the addendum at page [SI 
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In conclusion, in this paper we have suggested the possibility that the anomalies ob- 
served in the Gallium radioactive source experiments and the LSND experiment may 
be due to neutrino oscillations generated by the same large squared-mass difference 
^ m LSND+Ga- We have shown that, in the framework of the simplest 3+1 four-neutrino 

scheme that can accommodate also the ( v e — > ( z/^ T oscillations observed in solar and re- 
actor experiments and the ( z/^ —^v T oscillations observed in atmospheric and accelerator 
experiments, the short-baseline disappearances of v e and are due mainly to v e — > v s 

and ( v\ — > v s transitions, respectively. We have noticed that in the MiniBooNE exper- 
iment flavor transitions are effectively energy-independent and the disappearance of v e 
due to u e — ► v s transitions could affect the search for — > u e transitions, because of 
the natural v e contamination of the beam. Finally, we have remarked that the scenario 
under consideration could be tested in future experiments with pure v e and beams, as 
Beta-Beam and Neutrino Factory experiments. 



Addendum: First MiniBooNE Results 

After the completion of this paper, the MiniBooNE collaboration released their first re- 
sults concerning the search for — ► v e transitions generated by Attend [HI]- Since no 
significant excess of quasi-elastic charged-current v e events was observed above the cal- 
culated background for reconstructed neutrino energy E® E > 475 MeV, the two-neutrino 
Vn — > v e transitions generated by AmL SND are disfavored by the MiniBooNE data at 98% 
C.L. UZD- 

In the framework of the 3+1 four-neutrino scheme considered in this paper, the ab- 
sence of a signal due to — > v e appearance may be, at least partially, explained by a 
suppression of the background due to v e — > v s and — > v s transitions, as remarked after 
Eq. (I45p . In fact, the estimated number of v e events is 

N Ve = Nl + P n Nl + P v ^ Ve N v „ , (47) 

where iV^ and are, respectively, the estimated numbers of z/ e -induced and z/^-induced 
background events, and N v is the estimated number of N Ue in the case of full — > v e 
transmutation. In short-baseline experiments P Ve -+ Ve — 1 — P Ve -^, Us , as remarked after 
Eq. (1371) . and P„ ~ 1 — P v -> Us , as remarked after Eq. (T46l) . Moreover, the oscillation 
probabilities are practically constant in the MiniBooNE energy spectrum, as explained 
after Eq. (jlSj). 

From Table I of Ref . [HI] , adding the uncertainties in quadrature, we obtain 

= 229 ± 32.5 and JVj = 129 ± 17.0 , (48) 

for E^ E in the range 475 MeV < E® E < 1250 GeV. From the public information kindly 
given by the MiniBooNE collaboration on the Wetfl we obtain 

N Vll = 62851.2 + 250.7. (49) 

6 http : //www-boone . f nal . gov/f or_physicists/ april07datarelease/ 
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From Eqs. (l5]h (1551) and (j27j», we have 



P„ e ^ e = R = 0.88 ± 0.05 , (50) 

= l ~\ sin2 = °- 96 ± °- 02 , (51) 
Pu^u e = P^Z = (2.64 ± 0.81) x 10- 3 . (52) 



Then, we obtain 

*V^e < = 201 - 5 ± 30 - 8 , M ^ M = 123 ' 8 ± 16 ' 5 > *V-* = 165 - 9 ± 50 - 9 • 

(53) 

Comparing with iV^ in Eq. (|48|) . one can see that the estimated amount of z/ e -induced 
background is reduced by about 28 events as an effect of u e — > z/ s transitions. This 
reduction can compensate only partially the larger appearance signal due to — > v e 
transitions. 

The estimated and measured numbers of v e events are, respectively, 

N Ve = 491.3 ± 61.7 and jv^ iniBooNE = 380 ± 19.5 , (54) 

Hence, the 3+1 four-neutrino scheme considered in this paper is compatible with the 
results of the MiniBooNE experiment within 1.7 standard deviations. Although our 
scheme is clearly not favored by the MiniBooNE data, we think that further measurements 
are necessary in order to assess its viability. 
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